Distillable entanglement and area laws in spin and harmonic-oscillator systems 
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We address the presence of non-distillable (bound) entanglement in natural many-body systems. 
In particular, we consider standard harmonic and spin-1/2 chains, at thermal equilibrium and char- 
acterized by few interaction parameters. The existence of bound entanglement is addressed by 
calculating explicitly the negativity of entanglement for different partitions. This allows to individ- 
uate a range of temperatures for which no entanglement can be distilled by means of local operations, 
despite the system being globally entangled. We discuss how the appearance of bound entanglement 
can be linked to entanglement-area laws, typical of these systems. Various types of interactions and 
topologies are explored, showing that the presence of bound entanglement is an intrinsic feature 
of these systems. In the harmonic case, we analytically prove that thermal bound entanglement 
persists for systems composed by an arbitrary number of particles. Our results strongly suggest the 
existence of bound entangled states in the macroscopic limit also for spin-1/2 systems. 

PACS numbers: 03.67.Mn, 03.67.-a 



I. INTRODUCTION 

An ubiquitous scenario in Quantum Information Sci- 
ence consists of two (or more) parties sharing entangled 
quantum states and performing operations on them. Ide- 
ally the shared states are pure and maximally entangled. 
However, this is hardly the case in practice. As a mat- 
ter of fact, the available states are always mixed, due to 
unavoidable errors in the preparation stage or noise in 
the considered process. Entanglement is also very fragile 
in this sense, and this is certainly the main obstacle to 
actual applications of Quantum Information ideas. 

Entanglement distillation plays an important role to 
correct the degradation of entanglement in real situations 
[HQ- It consists on the application of local operations 
supplied by classical communication (LOCC) that allow 
the parties to extract maximally entangled states out of 
a bunch of mixed states. Unfortunately not all entangled 
states are distillable [3]. There are states p from which 
no LOCC strategy is able to extract pure-state entangle- 
ment, even if many copies of p are available. These states 
are known as bound entangled. 

The existence of bound entangled states was proven 
in Q by noting that any quantum state with a positive 
partial transposition (PPT) Q is non-distillable. Since 
there were already examples of PPT entangled states [j| , 
the connection between non-distillability and the positiv- 
ity of partial transposition automatically led to the first 
examples of bound entangled states. Later, strong evi- 
dence was provided for the existence of bound entangled 
states with non-positive partial transposition (NPPT) 
@. Independently of partial transposition, bound en- 
tangled states have a clear operational definition: an en- 
tangled state of n parties is bound entangled whenever 
the n parties cannot distill any pure-state entanglement 
out of it by LOCC. Although several examples of bound 
entangled states have appeared so far [7j, simple recipes 
to construct such kind of states are still lacking. It is thus 
not completely clear whether bound entanglement is ba- 



sically a mathematical artifact or it "naturally" appears 
in physically relevant quantum systems. The main mo- 
tivation of this work goes precisely along this direction, 
since we study whether bound entanglement is present 
in standard many-body quantum systems. By this, we 
mean systems i) characterized by a few interaction pa- 
rameters, ii) in thermal equilibrium with their environ- 
ment, Hi) consisting of a macroscopic number of particles. 

Very recently a few works have also addressed this 
question. In Ref. [1[ bound entanglement was detected 
in the thermal state of spin systems, consisting of up 
to 9 spins. More recently, three-qubit bound entangled 
states were obtained as the reduced state of the XY 
model in the thermodynamical limit 9]. In a previous 
contribution we have considered the problem of finding 
bound entangled thermal states in some specific system 
consisting of a macroscopic number of harmonic oscil- 
lators [lj|- The main goal of the present article is to 
give a detailed description of these results and extend 
them to other systems. Specifically, we directly calcu- 
late the entanglement in different partitions of thermal 
states for harmonic oscillators and spin systems and iden- 
tify a temperature range for which bound entanglement 
is present. In the harmonic case with nearest neighbor 
interaction, we first consider systems composed of hun- 
dreds of particles and then explicitly prove that bound 
entanglement indeed persists in the macroscopic limit, 
i.e. when n — » oo. Moreover we offer an explanation of 
the results in terms of the entanglement area law, a prop- 
erty satisfied by standard many-body systems [111 ]. Due 
to the general validity of area laws, we see that the pres- 
ence of bound entanglement is a common feature of these 
systems. We explicitly analyze a variety of different sys- 
tems, ranging from critical to non-critical situations and 
to spin chains characterized by a complex behavior of the 
entanglement in the ground state. In all these cases we 
see that there is a range of temperatures for which no en- 
tanglement can be distilled by means of LOCC, despite 
the system being globally entangled. 
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The paper is organized as follows. In the next Sec- 
tion we introduce the main intuition that elucidates the 
emergence of bound entanglement in thermal states, re- 
calling briefly some known results about entanglement- 
area laws. Then, in Sec. IIIII we present the analysis of 
harmonic systems. In particular, in Sec. IIIIBI numer- 
ical results regarding finite size systems are presented, 
whereas in Sec. MI CI we give the analytical calculations 
in the macroscopic limit for systems with nearest neigh- 
bor interactions. In Sec. HVl we show that similar results 
can also be obtained for spin systems, even if in this case, 
and due to numerical limitations, systems composed only 
by a small number of particles have been considered. We 
close the paper with some concluding remarks in Sec. |V] 



II. MULTIPARTITE BOUND ENTANGLEMENT 
AND AREA LAWS 

Let us present the main intuition behind our con- 
struction of bound entangled thermal states. Consider 
a quantum system of n particles described by a local 
Hamiltonian. For the sake of simplicity, in what follows 
we restrict our analysis to one-dimensional systems of n 
particles. A common property of these systems is that 
the ground-state entanglement obeys an area law. This 
means that the entanglement of a bipartite splitting of 
the system into two groups scales at most as the number 
of connections between the groups (i.e. the area which 
separates them). We recall that the entanglement for a 
bipartite splitting of pure states is uniquely quantified by 
the entropy of one of the subsystems [13], [lj] . In general, 
this quantity scales as the volume of the subsystem, and 
not as its boundary, for an arbitrary state belonging to 
the whole Hilbert space of the n particles. However, as 
said, this is not the case for the great majority of ground 
states of local Hamiltonian systems. 

Entanglement-area relations were first recognized in 
connection with the physics of black holes, for which the 
Von Neumann entropy scales as the surface at the event 
horizon [TEj . In recent years, the advances of entangle- 
ment theory allowed for the exploration of this behavior 
also for a variety of systems typical of condensed mat- 
ter physics (see Ref. [ll| and references therein). Re- 
markably, many analytical results have been found for 
harmonic systems for which concepts from the theory of 
Gaussian quantum states can be applied. As an example, 
a strict entanglement-area relation has been established 
for non-critical one-dimensional systems with finite range 
interaction, while logarithmic corrections appear in the 
critical case 0, [l7|, [ill ■ I n particular, the entanglement, 
as measured by the log-negativity [19j , has been proved to 
scale proportionally to the area. Similar behaviors have 
also been found for spin systems. Concerning thermal 
states, it has recently been demonstrated that the total 
amount of correlations (measured by the mutual infor- 
mation) in a bipartite split scale at most as the area [20| . 
This, in turn, gives an upper bound to the entanglement. 



Consider now a translationally invariant system com- 
posed by an even number n of particles (labeled from 1 to 
n) . Let us focus our attention on two different partitions 
of it, one in which a contiguous half of the particles be- 
longs to group A and the other half to B (we will refer to 
such kind of partition as half-half), and another partition 
in which the particles with even label belong to A and the 
rest to B (even-odd partition). Because of the area law, 
the entanglement saturates for sufficiently large n for the 
half-half partition, while it increases with n for the even- 
odd partition. In this configuration, it is reasonable to 
expect that, by increasing the temperature, the entangle- 
ment in the even-odd partition is more robust to thermal 
noise than in the half-half partition, and that this behav- 
ior is preserved for large systems. Denote by TA (T^°) 
the threshold temperatures at which the partial transpo- 
sition with respect to all half-half (even-odd) partitions 
becomes positive [2l|. Because of the area law, one can 
expect that Tjj£ h is strictly smaller than T^°. Thus, it 
emerges a range of temperatures for which the system is 
still entangled (as detected by the entanglement in the 
even-odd partition), nevertheless single particles cannot 
distill pure entanglement (as the half-half partitions be- 
come PPT). This is because, for any pair of particles, 
there is always a half-half partition for which they are 
in opposite sides and the partial transposition is posi- 
tive according to this splitting (remember that this is 
a sufficient condition for non-distillability [3]). In other 
words, bound entangled states are expected to appear un- 
der these conditions. The rest of the paper is devoted to 
verify and put on solid grounds this intuition for various 
relevant examples of many-body systems. 

Before proceeding, let us relate our finding to two re- 
lated works studying the presence of bound entanglement 
in systems with local interactions. First, in Ref. [8| the 
existence of thermal entangled states which are PPT with 
respect to any bipartition was proven for systems con- 
sisting of up to 9 spins. For such states, distillation of 
pure-state entanglement is impossible even if the parties 
arbitrarily combine together. As said above, we will in- 
stead show the existence of bound entangled states only 
with respect to fully local distillation procedures but for 
an arbitrary number of particles. Actually, for the har- 
monic systems considered here we explicitly prove that 
no entangled states PPT with respect to any bipartition 
exist. Second, in Ref. Q, the authors find bound entan- 
glement in the reduced three-qubit state of a macroscopic 
system. Here, we focus on the existence of bound entan- 
glement in the whole thermal state of the system. 



III. HARMONIC OSCILLATORS 

In this Section, we first introduce the harmonic sys- 
tems that we are going to consider in order to test the 
ideas exposed above. Then, in Sec. IIII B[ we perform 
numerical calculations for different types of interactions 
in systems of finite size. In Sec. IIII CI analytical results 
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are presented for the significant case of nearest-neighbor 
interaction in the macroscopic limit. 

A. The systems 

Consider a system composed of n harmonic oscillators, 
each one associated with position and momentum oper- 
ators Xi and pi respectively {i — 1, . . . , n), described by 
the Hamiltonian 

H = \ (e^+E^^) ■ w 

The diagonal elements of the matrix V describe the on- 
site interaction, while the non-diagonal terms give the 
coupling between oscillators i and j. Note that this 
Hamiltonian is quadratic in the canonical coordinates 
and the oscillators are coupled through their position de- 
grees of freedom. In this scenario both the ground and 
the thermal states turn out to be Gaussian, so they are 
completely described by their covariance matrix 7. Intro- 
ducing the vector S — . . . ,x n ,pi, . . . ,p n ), the latter 
is defined as follows 

lkl = Re (Tr{g[S k - S k ] [Si - §1}}) (2) 

where g denotes the density matrix of the state and 
Sk = Tr(gSk)- Considering the thermal state g = 
exp[-#/T]/Tr{exp[-i7/T]} at temperature T, the cor- 
responding covariance matrix is given by 22]: 

T (T) = [V-^ 2 W(T)} 8 [V l ' 2 W{T)] , (3) 

where 

W(T) = t n + 2[exp{V 1 ' 2 /T)-t n ]- 1 (4) 

and 1„ denotes the n x n identity matrix. In the ground- 
state case W(0) is given by the identity matrix and thus 

7 (o) = v- 1 / 2 V 1 ' 2 . 

The entanglement properties of the ground and ther- 
mal states corresponding to Hamiltonian UJ were first 
studied by Audenaert et al. |22l| . There, an analytical 
expression for the entanglement (quantified by the log- 
negativity [l9[ ) between two complementary groups of 
oscillators, A and B, was given in terms of the covari- 
ance matrix of the state, which can be written, in turn, 
only in terms of the matrix V. Then one gets the gen- 
eral formula for the log-negativity of a thermal state at 
temperature T: 

17-1 

Ei =^log 2 {max[l,A fc (Q)]}, (5) 

fc=0 

where Q = P uj~ Puo+ and w ± = W^T)" 1 !^. We 
denoted by {Afe[Q]}J?~Q the spectrum of the matrix Q, 
whereas P is an n x n diagonal matrix with the i-th en- 
try given by 1 or —1 depending on which group, A or 



B, oscillator i belongs to. This study was later extended 
in Ref. [l6[ where an area law for ground-state entangle- 
ment was proven. As far as for thermal states, an upper 
bound for the entanglement in terms of the number of 
connecting points in a given bipartition was also estab- 
lished |18| , but the lack of a lower bound limits us to get 
a strict area law in this case. 

In the case of a harmonic chain with nearest-neighbor 
interactions and periodic boundary conditions, the sys- 
tem is described by the Hamiltonian (UJ) with a circulant 
potential matrix V given by 

V n =circ(l,-c,0,...,0,-c). (6) 

The system is defined for < c < 1/2 and it is equiv- 
alent to a chain of harmonic oscillators coupled with a 
spring-like interaction. A variety of physical systems can 
be modeled by such interaction, going from vibrational 
degrees of freedom in crystal lattices and ion traps to free 
scalar Klein-Gordon field. We recall that its ground-state 
ent ang lement exhibits a critical behavior when c — > 1/2 

[2, Hall. 

In order to test the general validity of the connection 
between bound entanglement and area laws mentioned 
in Sec. [Til we have considered also other systems. As 
an example we will report here some results concerning 
the next-to-nearest interaction defined by the potential 
matrix 

V nn = circ(2 + 4/i 2 , -An, 1, 0, . . . , 0, 1, -Afx) . (7) 

This system has been shown to be critical (gapless) and 
to violate the area law for < fi < 1 in the ground state. 
Nevertheless, as already mentioned, no violation of the 
area law is possible for non-zero temperatures in view of 
the results in Ref. [20| . 

B. Numerical results 

Let us first consider the nearest-neighbor interaction 
given by Eq. ©. We used Eq. to compute the log- 
negativity for the even-odd and the half-half partitions 
\2l\ for different temperatures and number of particles. 
As already recalled, in Ref. [22j it is shown that a strict 
area law holds for this system in the ground state. Ac- 
tually, our calculations show that the log-negativity fol- 
lows a strict area law for non-zero temperatures as well. 
We depicted in Fig. [TJ the log-negativity as a function 
of the number of particles, for fixed coupling and differ- 
ent temperatures. One can clearly see that, apart from 
a transient for small n, the log-negativity increases lin- 
early with n for the even-odd case, while it saturates 
for the half-half partition. Furthermore, the change of 
the system temperature just affects the rate at which en- 
tanglement increases with n for the even-odd partition 
and the entanglement saturation value for the half-half 
partition. The validity of a strict area law has then a 
remarkable consequence: being the log-negativity depen- 
dence on n the same (linear or constant) for all T, then 
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FIG. 1: Log-negativity Ei as a function of the number of 
oscillators n for half/half (left panel) and the even/odd (right 
panel) partitions. Temperatures are T = 0.45 (solid line), 
T = 0.4 (dashed line), and T = 0.35 (dotted line), and c = 
0.4 (the same behavior is found for other values of c). For 
T = 0.45 the state is bound entangled. 



Ei goes to zero as T increases independently of n. In 
other words, the threshold temperatures TA and 
should be independent of the size of the system (apart 
for small n, for which the area law shows a transient). 
Regarding the presence of bound entanglement, we see 
in Fig. [T] that, e.g. , for T = 0.45 the log-negativity in 
the half-half partition is zero, meaning that single parti- 
cles cannot distill, nevertheless the system is entangled, 
as shown by the non-zero log-negativity in the even-odd 
partition. For the reasons exposed in Sec. |TT] we then 
conclude that for that temperature the system is bound 
entangled. 

An exhaustive evidence of the presence of bound en- 
tanglement for any coupling constant is given in the T—c 
diagram of Fig. O There we see that the threshold 
is strictly larger than TRf 1 for any c, indicating the exis- 
tence of a bound entanglement region as soon as c =/= 0. 
The calculations are in this case performed for systems 
composed by 800 oscillators. The range of temperatures 
T$ h < T < for which bound entanglement is guar- 
anteed is seen to increase when the system approaches 
the critical point c = 0.5. 

As said, the fact that all the half-half partitions are 
PPT implies that no entanglement can be distilled by 
means of local operations performed on each particle. 
However, it does not exclude the possibility that each os- 
cillator is entangled with the others [25| . Interestingly, we 
also found a region of temperatures for which none of the 
oscillators is entangled with the rest of the chain, yet the 
global state is entangled (shaded region in Fig. ^ . This 
is shown by calculating the log-negativity for the parti- 
tions 1 : n — 1, i.e. one particle versus the others. Since 
the partition consists of 1 versus n — 1 modes, the PPT 
condition turns out to be sufficient for separability [26| . 
We have calculated also the threshold temperatures cor- 
responding to a generic contiguous splitting §— to : §+to 
(with < m < 5 — 1). It can be seen that these threshold 
temperatures increase with m, that is, the entanglement 
in the bipartition 1 : n — 1 is the more robust to thermal 
noise among all these splittings. 

In order to analyze the dependence of the results above 
on the system size, we have plotted the threshold temper- 
atures as a function of n (see the inset of Fig. In ac- 
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FIG. 2: Threshold temperatures above which the log- 
negativity is zero in the even-odd (T^°, solid line) and half- 
half (T t h , dotted line) partitions. The thresholds are plotted 
as a function of the coupling constant c for the Harmonic 
chain with nearest-neighbor interactions composed by 800 os- 
cillators. The intermediate line (dashed) shows the temper- 
ature T^ n ~ above which the state is separable in all the 
partitions 1 : n - 1. Inset: (solid line), T& h (dotted 

line), and T^™" 1 (dashed line) as a function of the number 
n of oscillators composing the system (log-lin scale). The os- 
cillators interact via nearest-neighbors couplings with c = 0.3 
(the same behavior is found for other values of c). The gap 
Tth° — ?th i s seen to remain constant with the size of the 
system, apart from an initial transient. 

cordance to what was discussed, we see that the threshold 
temperatures remain constant. As a consequence, also 
the gap T t e jj° — XA , which determines the temperature 
range in which bound entanglement is present, is inde- 
pendent of the size of the system. These results strongly 
suggest that bound entanglement can also be observed 
in the macroscopic limit, an issue on which we will come 
back in the next session. 

Due to the generality of the area law, the results ob- 
tained for the nearest-neighbor model above are expected 
to be valid in a variety of different cases. As an example 
we report here the next-to-nearest interaction V nn intro- 
duced in the previous subsection. In the T — \i diagram 
of Fig. [3] we can see that there is again a wide range of 
temperatures for which bound entanglement is present, 
as detected by the coexistence of nonzero log-negativity 
in the even/odd partition and zero log- negativity in the 
half/half one. Notice that the absence of an energy gap 
for fi < 1 does not affect qualitatively the presence of 
bound entanglement [27|- As said above, this can be re- 
lated to the fact that an area law also holds for gapless 
systems in case of non-zero temperatures. 



C. Analytical results 

As pointed out in the previous section, the validity of 
a strict area law suggests both the independence of the 
threshold temperatures from the size of the system as well 
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FIG. 3: T&° (solid line), T^' 1 " 1 (dashed line), and T^ h (dot- 
ted line) as a function of the coupling constant /i for the Har- 
monic chain with next-nearest-neighbor interactions in Eq. ([7} 
composed by 200 oscillators. 

as the preservation of thermal bound entanglement for 
macroscopic systems. In this section we will analytically 
prove that this is actually the case. In particular, we 
will establish, for the paradigmatic system in Eq. ([5]), a 
formula for Ei in the even-odd partition and an upper 
bound for in the half-half partition. Both will be 
obtained in the macroscopic limit, i.e. for n — > oo. The 
remaining part of this Section is quite technical. The 
reader not interested in the details can skip the proof, 
whose results are for convenience summarized in Fig. [?] 
Let us start with the case of the even-odd partition, for 
which we can analytically prove the linear increase of the 
log-negativity with the system size, as well as obtain the 
threshold temperature as a function of the coupling 
c. The proof involves the calculation of the spectrum of Q 
in Eq. ([5]), and it follows closely the one given in Ref. [22] 
for the ground-state case. Consider the Hamiltonian ([1]) 
and a generic potential with translational invariance sym- 
metry and periodic boundary conditions. Then, V is a 
circulant matrix and the matrices can be diagonalized 
by a discrete Fourier transformation implemented by 

Sl kt l=expikl J /Vn, (8) 

with < fc, I < n—1. Namely, one has that f2w ± f2^ = D ± 
with D± = Skidf and df = A^ 1/2 taah(\/A~k/2T) where 
Afc are the eigenvalues of V. Concerning the matrix P, 
one has that: 

nM, - p -(^,V)- < 9 > 

The spectrum of Q coincides then with the one of 
PD~ PD + , which in turn is straightforward to calcu- 
late due to its block diagonal structure. The eigenval- 
ues of Q are hence given by d^d^ + , 2 and d^ +n ^ 2 d^ , 
for k = 0, . . . , n/2 — 1. Specializing these results to the 
nearest- neighbor case in Eq. ©, one has that the eigen- 



values of V n are given by Afc = 1 — 2c cos (2irk/n). Defin- 
ing the function 

(10) 

one has that the eigenvalues of Q that can contribute to 
the log-negativity (i.e. , that can be larger than one in 
dependence of the parameters) are given by /(fc, n, c, T), 
with double multiplicity and fc = Q, . . . , n/4 (for n multi- 
ple of 4, n > 4). The function /(fc, n, c, T) is non increas- 
ing in fc for the range of interest (0 < c < 1/2 and T > 0), 
and in particular it attains its maximum for fc = 0. As a 
consequence, the log-negativity of the even-odd partition 
is different from zero when the temperature T is such 
that /(0, n, c, T) > 1, namely: 

In particular, the curve /(0, n, c, T) = 1 gives the thresh- 
old temperature T^°. It coincides with the one depicted 
in Fig. [5] (solid line) and it is worth noticing that it does 
not depend on the total number of particles n (see also 
the solid line in the inset of Fig. [2]), i.e. it holds also 
in the macroscopic limit. For temperatures below T t e jj° 
there exists a k(n,c,T) such that f(k,n,c,T) > 1 for 
fc < k(n,c,T), which in turn gives rise to the following 
expression for the log-negativity: 

k(n,c,T) 

Ei= log 2 {/(fc, n,c,T)}. (12) 

Following again Ref. [22j , we have that for large n the sum 
over fc can be replaced by an integral over x = 2irk/n : 

n r {c ' T) 

Ei = ^J dx \og 2 { f(x,c,T)}. (13) 

Since f(k,n,c,T) in Eq. (fTO)) actually depends only on 
the parameter x, we have that both the integrand and 
the upper bound x(c, T) turn out to be independent of 
n. As a consequence the log-negativity grows linearly 
with the system size in the macroscopic limit also for 
non-zero temperatures. 

Regarding the half-half partition, we were not able to 
find an analytic expression for the log-negativity in the 
macroscopic limit. Nevertheless it is possible to find an 
exact upper bound for the threshold temperature T^ h . 
As said, this allows to identify an interval of temperatures 
for which the presence of bound entanglement can be 
guaranteed also in the macroscopic limit. We proceed as 
follows. Let us define the matrix 

n/2-l n -l 

Xij=U^ 22 X! (Sikfijh + SjkSih) (14) 
fc=0 h=n/2 



6 



where Sij denotes the Kronecker delta. According to 
Ref. [TH, the log- negativity in the half-half partition is 
zero when the following inequality is satisfied: 

[W(T)] 2 + 2 maxi|Ai[Xijj + ]| < 1 , (15) 

where the minimum eigenvalue of W(T) is given by: 

/T+2S/T + 1 



Xmin[W(T)] = 



(16) 



Recognizing that the second term in the left hand side 
of Ineq. (fT5|) is twice the spectral radius r{Xu> + ) of the 
matrix Xlo + , we can use any matrix norm to bound it 
from above [28|. In particular, let us consider the maxi- 
mum row sum matrix norm || • H^, defined for a generic 
matrix A as: 



max, 



3 



(17) 



(in what follows, we will omit the subscript for the sake 
of convenience). An upper bound for r(Xu> + ) is then 
given by r(Xu> + ) < \\Xoj + \\, from which it follows that: 

r(Xcu+) <\\X\\\\u + \\ . (18) 

Before proceeding in bounding \\X\\ and ||w + ||, let us 
write explicitly their elements in the macroscopic limit. 
Recall that ui^ are circulant matrices, hence completely 



specified by their first row u> 



± _ 



circ(i; , . . . ,v n _i_), 



whose elements are in turn given by the discrete Fourier 
transform of the eigenvalues d k introduced above: 



E 

fe=0 



exp 



2ni. 



-kl 



(19) 



In the limit n — > oo the discrete sum above yields the 
Riemann sum for the following integral in the argument 
x = lixkjn: 



1 

2^ 



dx d ± (x)e 



ixl 



(20) 



that is, u, are the Fourier coefficients of the periodic 
functions d ± (x): 



d ± {x) = (1 - 2ccosa;) ±1/2 tanh( 



VI — 2ccosx 
2T 



)• (21) 



Notice that d (x) are smooth functions, which implies 
an exponential decay of \v^\ as a function of I. In fact, 
for any integer s, by integrating by parts s times we have: 



< 



1 

2^ 
1 I 
2ttF 



1 



(ay 



2tt 



dx e lxl ■ ' 



P 1 " 



(22) 



where Ct~ 



Jo* dx 



dx 



Let us now bound first 



Being u> + a circulant matrix, there is no need to 



look for the maximum in (|17p , as the rows are composed 
by the same elements. Then one can write, for any integer 
to, ||w + || = 5+ + £+, where we defined the partial sum 
and the residual term, respectively, as follows: 



'3 n 



E i 

=m+l 



(23) 



In order to obtain a bound on ||w + || one can now fix to, 
calculate explicitly and bound £^ from above. This 
latter step can be achieved by exploiting again the prop- 
erty of d + (x) and integrating by parts s times, leading 
to: 



e+ < ^C(«,m + 1), 



(24) 



where £(s, to + 1) is the generalized Riemann zeta func- 
tion. Summarizing, one has that for any integer s and 
m 



k + ||<S+ + ^C(s,m + l). 

7T 



(25) 



the bound becoming tighter for large to and s. 

Let us now determine an upper bound for \ \X\\. First, 
for any finite n, the matrix X has a bipartite symmetric 
structure of the following form: 



X = 



B 
B 



where B is a 5 x ^ symmetric Toeplitz matrix. Then \ \X\\ 
coincides with ||_B|| and in order to evaluate the latter 
we need to determine which row of B gives rise to the 
maximum in Eq. (|17p . Fortunately, in the macroscopic 
limit n — * oo the periodic boundary conditions can be 
disregarded and the last row of B determines its norm. 
This is because going from the last to the next to the 
last row we simply remove the term from the sum 
in Eq. (J7J), and so on for the other rows. It can be seen 
that this is the case also taking rigorously into account 
the periodic boundary conditions. By writing explicitly 
matrix B: 



t 



B 



i/2 



a/2-1 



2/2-1 



V 



i/2-l 



n/2 
V~ 



/2 



/ 



and exploiting its symmetries, we see that our consider- 
ations can be restricted to the last n/4 rows (or, equiva- 
lently, to the first n/4 rows). Define for a generic row the 
following sum (J — 1, . . . , n/4 denotes the rows in reverse 
order, from the last to the j-th one): 



R; 



n/2-] 

E 



n/2 

E 



(26) 



\+l-3 
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such that \\B\\ = maxj Rj. The difference between Rj 
for a generic row and i?i can be bound from above as 
follows: 

re/2 j n/2 

Ri-Ri= E ITI-EITI< E (27) 

which implies that 

re/2 

maxj(flj -i?i) < E l w fl- ( 28 ) 
i=? +1 

Since the sum above goes to zero for n — > oo [as can be 
explicitly seen by exploiting the exponential decay of \vf \ 
in Eq. (|22[) ] . we have that the error made by considering 
the last row of B as the one that gives the maximum in 
evaluating ||B|| can be made as small as one may wish. 
In other words, we can identify \\X\ \ = R\. 

Now, we can proceed as we did for | \to + \ |. In particular 
defining again the partial sum as 

m 

S- = ElVh (29) 
1=1 

we obtain that for any integer s and m 

\\X\\<S- + ^C(s,m+l), (30) 

where again S 1 " can be calculated explicitly and the 
bound becomes tighter for large m and s. 

Summarizing, considering Eqs. (|15p . (flB)) . (I25)) , and 
pop , we have shown that, in the macroscopic limit, the 
log-negativity in the half-half partition is zero when the 
following inequality is satisfied: 

2 (s+ + ^C(s, m + 1)) (s~ + ^((s, m + 1)) 

+ ( e^ + l j ° (3 " 

Based on the formula above and on Ineq. (dJ), we de- 
picted in Fig. 3] the region in the c — T plane for which 
bound entanglement is present in the macroscopic limit 
(shaded region). We see that for any coupling c we can 
guarantee that there is an interval of temperatures for 
which the negativity in the half-half partition is zero, 
nevertheless the state is entangled. 

Finally, let us mention an issue concerning the proce- 
dure used in order to derive Ineq. We considered 
a finite system composed by n particles and calculated 
explicitly T^°. We have seen that it does not depend on 
n and so argued that it has to be valid also in the limit 
n — ► oo. On the other hand, one could have considered 
first the limit n — > oo and then analyzed a distinguished 
region of the system for increasing size [as we did in de- 
riving Ineq. (|3ip ]. The two procedures are equivalent, 
since no convergence problem appears in the spectrum of 
V n for n — > oo. 



0.6F 




0.1 0.2 0.3 0.4 0.5 

c 



FIG. 4: The solid line represents in the macroscopic limit, 
accordingly to Ineq. {TTJ : below this line the state of the sys- 
tem is entangled. The dashed line represents the upper bound 
to T^ h as given by Ineq. (|3ip (e.g., we considered m = 10 and 
s — 3): above this line the log- negativity in the half-half par- 
tition is zero. In the shaded region then we can guarantee the 
presence of bound entanglement in the macroscopic limit. 

IV. SPIN SYSTEMS 

As already mentioned, we go beyond the harmonic os- 
cillator case and also investigate the existence of bound 
entangled thermal states in spin systems. Consider the 
Hamiltonian 

n n 

Hxx = -JE«<+i + ff M+i) + B E<> ( 32 ) 

i— 1 i 

that is, n spin-one-half particles with nearest-neighbor 
XX interactions with coupling constant J and subject 
to a transverse magnetic field B. For this system we also 
calculate the thermal-state entanglement, measured now 
by the negativity En, for different partitions. All the ob- 
tained results are consistent with the previous reasoning: 
there is a temperature range for which the negativity in 
the half-half partitions is zero, nevertheless the system is 
still entangled as proven by the negativity in the even-odd 
partition. Although in this case we are not able to deal 
with big systems, our results suggest an entanglement 
area law for non-zero temperature and the independence 
of the gap T t e jj° — T^ h with the system size, as found in 
the harmonic oscillator case. 

In Fig. [5] we plot the negativity for the even-odd (full 
symbols) and for the half-half (empty symbols) parti- 
tions, for two different values of the temperature. For a 
constant number of bonds between the partitions as in 
the case of the half-half splitting the negativity does not 
depend on the system size n, whereas for the even-odd 
partition it increases with n. This behaviour resembles 
the previous study for harmonic oscillators (see Fig. [1]), 
and the same reasonings apply here to show the presence 
of bound entanglement. 

In order to study the PPT condition in a system with 
n = 10 particles for different values of the coupling con- 
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FIG. 5: Negativity in the even-odd (full symbols) and half- 
half (empty symbols) partitions for the thermal states of the 
Hamiltonian (|32|l with B — 1.9 and J = 1. The temperatures 
for each partition are set to be T — 2 and 2.6 (see legend). 
In the even-odd partition we can clearly see an increase of 
the negativity with respect to the system size, whereas it sat- 
urates in the half-half partition. For T = 2.6 we see that 
the state is bound entangled (the negativity in the half-half 
partition being zero). 




FIG. 6: Threshold temperatures above which the negativity 
is zero in the even- odd (Tg°, solid line), 1 : n - 1 (T t ^ n_1 , 
dotted line) and half-half (T t ^ , dashed line) partitions. We 
plot the threshold temperature as a function of the coupling 
parameter J of the Hamiltonian (|32[l with n = 10 and B = 
1.9. Inset: Temperature above which the negativity in the 
even-odd (T£°, solid line), 1 : n- 1 (T^ n_1 , dotted line) and 
half-half (T^ , dashed line) partitions is zero as a function of 
the number n of spins with Hamiltonian (|32|l with J = 1 and 
B = 1.9. 



stant J, we plot the corresponding separability temper- 
ature in Fig. [B] for three different bipartite splittings. 
In all the range of values of J the separability temper- 
ature corresponding to the even-odd partition TS° is 
strictly higher than T t ™~ and T^ h , corresponding to 
the 1 : n — 1 and half-half partitions respectively. The 
corresponding threshold temperatures for PPT are shown 
in the inset of Fig. [S] for increasing system sizes. For 
n > 6 the gap between these temperatures appears con- 
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FIG. 7: Negativity in the even-odd (solid line) and half-half 
(dahed line) partitions for the thermal states of Hamiltonian 
l|32l) with n = 10, J = 1 and T = 0.1 as a function of B. For 
B > 2 the ground state is the separable state 00 . . . 0), but 
some amount of entanglement still appears due to the mixture 
with entangled excited states. 



stant (we only compute these values for systems with an 
even number of particles for symmetry reasons). 

We can tune the local field B of the hamiltonian 
Eq. ((321) in such a way that the the system ground state 
exhibits some interesting entanglement properties. For a 
local field B = 1.9 and J = 1 the ground state of the 
system is a generalized W-state of n spins [2£| , whereas 
for B > 2 the ground state is the product state |00 . . . 0). 
In Fig. [7] we plot, for different values of the field B and 
T = 0.1, the negativity for the even-odd and half-half 
partitions. Notice the presence of the characteristic dips 
due to the mixing of different entangled ground states 
[30| , In particular, we can see that for B < 2 the even- 
odd partition has a higher degree of entanglement than 
the half-half partition. As a consequence, the main points 
of our intuition for the emergence of bound entanglement 
are still valid, also in the presence of complex behaviors 
like the ones shown in Fig. [JJ 

For low temperatures and values of the local field 
B > 2 the ground state (separable) dominates the sys- 
tem and therefore the thermal state is separable for all 
the partitions. For higher temperatures the state turns 
out to be a combination of the separable ground state and 
some excited states that can be entangled [3l|. In this 
condition entanglement can be increased by raising the 
temperature. This is shown in Fig. [51 where we plot the 
negativity for fixed field B = 2.3 as a function of the tem- 
perature for the even-odd and half-half partitions. We 
can see again that the negativity in the even-odd parti- 
tion is higher than the respective values for the half-half 
partition. In particular, in the range 2.5 < T < 3.2 
the half-half partition is PPT while the negativity in 
the even-odd partition is non-null. Thus, also in this 
quite peculiar conditions - in which the ground state has 
-Ejv = - we see, for increasing temperatures, that i) 
when the system becomes entangled (i.e. , for small T 
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FIG. 8: Negativity in the even-odd (solid line) and half-half 
(dahed line) partitions for n = 10, J = 1 and B = 2.3 for 
spin-one-half particles with Hamiltonian (|32|l . As expected 
the negativity in the even-odd partition is greater than the 
corresponding to the half-half partition, and bound entangle- 
ment appears in the temperature range 2.5 < T < 3.2. 

greater than zero) the negativity is higher in the even- 
odd partition and ii) before the system becomes separa- 
ble (i.e. , for higher T) it passes through a non-distillable 
region. This behavior again supports the general validity 
of our considerations. 

Finally, we mention that we have found an analogous 
behavior in systems subjected to the Heisenberg interac- 
tion 

n n 

H X xx =E(«+i+^ CT m+< <7 ?+i)+ s E^ ( 33 ) 

i—1 i 

for which a study of the entanglement properties of two 
neighboring spins was shown in 'Mi . In this model a 
value of the local field B > 4 produces the ground state 
1 00 . . . 0). From this point the results and considerations 
for this model are essentially equivalent to those corre- 
sponding to the XX model. 

V. CONCLUDING REMARKS 

In closing this paper, we point out a few remarks re- 
garding the obtained results. First, we have shown the 
existence of thermal bound entanglement by considering 
local (n-party) distillation strategies. One can also pose 
the question whether fully PPT entanglement can ex- 
ist in these systems, that is entangled states such that 
all bipartite partitions arc PPT. Actually, as shown in 
Ref. Q, fully PPT entanglement can be found for spin 
systems, at least up to 9 qubits. On the other hand, our 
results, when combined with those of Refs. [H,[33[, prove 
that fully PPT entangled states cannot be obtained for 
the harmonic systems studied here. In Ref. [32| bounds 



are derived for the threshold temperatures for full sepa- 
rability of various harmonic systems. Remarkably, these 
upper bounds are tight in case of systems with sufficient 
translational symmetries, as in the case of the nearest- 
neighbor interaction in Eq. ([6]). It is then possible to 
depict a kind of phase diagram in the T — c plane sep- 
arating two regions, one in which the state is separable 
and the other in which quantum correlations are present 
in some form. Such a threshold turns out to coincide 
with the one given by Ineq. for the zero negativity 
in the even/odd partition. Thus, the thermal states of 
this system become PPT and fully separable for the same 
temperature. We mention here that such a result holds 
also for other systems, for example for the one given in 
Eq. ([7]). However, in these cases we had to approach the 
problem numerically (based on the results of Ref. (33|) 
and consider only finite systems. 

Finally, let us stress that our results enrich the above 
mentioned phase diagram of Ref. [32j , in the sense that 
we can distinguish a third region in it where quantum 
correlations are present but they are not distillable. In 
particular, referring to Figs. [2] and [4] we have that, by 
rising the temperature at fixed coupling, the system goes 
from a genuinely entangled state (in which all the bipar- 
titions are NPPT) to a bound entangled state and finally 
becomes just classically correlated. As already discussed, 
we have found that such a behavior is common to all the 
harmonic systems we took in consideration. For spin 
systems, a similar phase diagram holds, at least for the 
finite-size systems studied here. 

In conclusion, we have shown the existence of thermal- 
state bound entanglement for several quantum many- 
body systems. We explicitly considered systems com- 
posed by hundreds of oscillators and observed in gen- 
eral an area law for thermal states. This findings turned 
out to be crucial in proving the existence of bound en- 
tanglement also in the macroscopic limit for harmonic 
systems. As a consequence, and by virtue of the gener- 
ality of area laws, these results have been shown to be 
valid for a variety of other systems, as well as in the pres- 
ence of complex behaviors of the entanglement. All these 
considerations strongly support the existence of thermal- 
state bound entanglement in the macroscopic limit also 
for finite-dimensional systems. 
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